Fast articulated motion tracking

ABSTRACT

The present technology relates to a computer-implemented method for tracking an object in a sequence of multi-view input video images comprising the steps of acquiring a model of the object, tracking the object in the multi-view input video image sequence, and using the model.

RELATED APPLICATIONS

This application is a continuation of, and claims priority to, International Application No. PCT/EP2012/055650 (International Publication Number WO 2012/156141), entitled “Fast Articulated Motion Tracking,” and having an International filing date of Mar. 29, 2012. International Application No. PCT/EP2012/055650 claims priority benefits from European Patent Application No. 11166259.9, filed May 16, 2011, and U.S. Provisional Patent Application No. 61/486,380, filed May 16, 2011. International Application No. PCT/EP2012/055650, European Patent Application No. 111166259.9, and United States Provisional Patent Application No. 61/486,380 are hereby incorporated by reference herein in their entireties.

BACKGROUND

The present invention relates to fast articulated motion tracking.

One of the fundamental problems in computer vision is estimating the 3D motion of humans. Motion capture is an essential part in a wide range of modern industries, ranging from sport science, over biomechanics to animation for games and movies. The state-of-the-art in industrial applications is still marker-based optical capture systems, which enable accurate capture at the cost of requiring a complex setup of cameras and markers. Marker-less methods that are able to track the motion of characters without interfering with the scene geometry, rely on pose detection based on some image features. These approaches, however, assume that the poses have been previously observed during training.

Most related to the present invention is a real-time tracking system called Pfinder proposed by Wren et al. [C. Wren, A. Azarbayejani, T. Darrell, and A. Pentland. Pfinder: Real-time tracking of the human body. TPAMI, 19:780-785, 1997]. It models the human by 2D Gaussians in the image domain and represents the appearance of each blob by an uncorrelated Gaussian in the color space. The background is modeled by Gaussians in the color space for each image pixel. Pose estimation is finally formulated as 2D blob detection, i.e., each image pixel is assigned to the background or to one of the human blobs. The final 2D pose is obtained by iterative morphological growing operations and 2D Markov priors. The approach has been extended to the multi-view case in [S. Yonemoto, D. Arita, and R. Taniguchi. Real-time human motion analysis and ik-based human figure control. In Workshop on Human Motion, pages 149-154, 2000] where the blobs are detected in each image and the 3D position of the blobs is then reconstructed using inverse kinematics.

Other approaches to Human pose estimation without silhouette information combine segmentation with a shape prior and pose estimation. Graph-cut segmentation has been used as well as level set segmentation together with motion features or an analysis-by-synthesis approach. Handheld video cameras and a structure-from-motion approach are also used to calibrate the moving cameras. While these approaches iterate over segmentation and pose estimation, the energy functional commonly used for level-set segmentation can be directly integrated in the pose estimation scheme to speed-up the computation. The approach, however, does not achieve real-time performance and requires 15 seconds per frame for a multi-view sequence recorded with 4 cameras at resolution of 656×490 pixels.

Implicit surfaces have been used for 3D surface reconstruction. More particularly, a human model comprising a skeleton, implicit surfaces to simulate muscles and fat tissue, and a polygonal surface for the skin have been used for multi-view shape reconstruction from dynamic 3D point clouds and silhouette data. The implicit surfaces are modeled by Gaussians. Since the estimation of the shape and pose parameters is performed by several processes, the whole approach is very time consuming and not suitable for real-time application. To increase the reliability, motion priors can be used to improve model fitting. An implicit surface model of a human has been matched to 3D points with known normals.

It is therefore an object of the invention, to provide a method for acquiring a model and for tracking the pose of an actor in a digital video that is fast and efficient.

SUMMARY

This object is achieved by the methods and a device according to the independent claims. Advantageous embodiments are defined in the dependent claims.

According to the invention, a human model may be represented by a set of spatial Gaussians instead of making use of an explicit surface modeled by geometric primitives or a detailed triangle mesh. To this end, a person-specific model may be created from a sparse set of images. The model may comprise a kinematic skeleton that defines the degrees of freedom (DoF) of the human model and a statistical model that represents the shape and appearance of the human. In order to achieve real-time performance for pose estimation from multiple views, Sums of 3D Gaussians (SoG) may be used for the statistical model. Furthermore, the image may also be represented as a sum of spatial 2D Gaussians that cover image blobs that are consistent in color. Based on the representation of the model and the images as Sums of Gaussians, the invention introduces a novel formulation of the model-to-image similarity, whose analytical solution may be solved very efficiently with a gradient ascent algorithm.

More particularly, it allows a nearly real-time approach that does not rely on extracted silhouettes or training data.

Also, the inventive approach does not rely on a statistical background model to detect 2D blobs in the image, but uses 3D spatial Gaussians for the human model and 2D spatial Gaussians for the images to introduce an objective function for the model-to-image similarity that can be very efficiently maximized.

Finally, the approach does not rely on 3D data clouds or silhouette data for tracking, but models the similarity between the model and the unsegmented image data directly.

In some embodiments, the present technology provides a computer-implemented method for tracking an object in a sequence of multi-view input video images. The method can comprise the steps of acquiring a model of the object, and tracking the object in the multi-view input video image sequence, using the model.

Some embodiments present a system for tracking an object comprising a range of multi-view video cameras, a memory for storing a model of the object, and a processor for tracking the object based on the model.

In certain aspects, the present technology provides a computer-implemented method for acquiring a model of an object. The method can comprise the steps of acquiring a plurality of images of example poses of the object, estimating the model based on the acquired images, and storing the model in a computer-readable medium. The present technology also provides a computer-readable medium comprising a model, wherein the model is acquired by performing the aforementioned method.

BRIEF DESCRIPTION OF THE DRAWINGS

The patent or application file contains at least one drawing executed in color. Copies of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.

FIG. 1 shows an overview of a method according to an embodiment of the invention.

FIG. 2 shows, from left to right, a default skeleton, a default body model and an actor-specific model according to an embodiment of the invention.

FIG. 3 is a schematic diagram illustrating self-occlusion handling according to an embodiment of the invention.

FIG. 4 shows an example of SoG image approximation according to an embodiment of the invention. Left: Input image. Right: Quad-tree structure with average colors used to generate the SoG. Each square is represented by a single Gaussian.

FIG. 5 shows an example of estimating an actor specific model from example pose images according to an embodiment of the invention. Left: Single segmented input image of the multi-view sets for each pose. Right: Resulting actor-specific body model after optimization and color estimation.

FIG. 6 shows tracking results of a method according to an embodiment of the invention as skeleton overlay over 12 input images from different camera views.

DETAILED DESCRIPTION

FIG. 1 shows an overview of a method according to an embodiment of the invention.

In a pre-processing step 110, a low number of manually segmented multi-view images showing example poses may be used to estimate an actor specific body model. The human body model comprises a kinematic skeleton and an attached body approximation.

The acquired model may then be used for tracking the articulated motion of the actor from multi-view input videos. More particularly, the performance of an actor may be captured with n_(cam) synchronized and calibrated video cameras. Each input image acquired in step 115 may be converted into a body model representation in step 120.

Tracking may start with the estimated pose of the model in the previous frame, and may optimize the parameters such that the overlap similarity between the model and image SoG at the current frame is essentially maximized.

FIG. 2 shows, from left to right, a kinematic skeleton model, a default body model and an actor-specific model according to an embodiment of the invention.

The kinematic skeleton model according to the present embodiment of the invention comprises 58 joints, modeling a detailed spine and clavicles. Each joint may be defined by an offset to its parent joint and a rotation represented in axis-angle form. It may also feature an allowable joint limit range l_(l) to l_(h). The model comprises a total of 61 parameters Λ, 58 rotational and an additional 3 translational. The kinematic skeleton model further comprises a separate degree of freedom (DoF) hierarchy, having nDoF pose parameters Θ. In the present embodiment of the invention, a DoF hierarchy comprising n_(DoF)=43 pose parameters was used.

A relationship between Θ and Λ may be modeled as a linear mapping given by an 61x_(nDoF) matrix M:

Λ=MΘ  (1)

where each entry of M defines the influence weight the parameters of Θ have on the joint angles Λ. This construction allows the kinematic skeleton model to reproduce a natural deformation of the spine, as a single DoF can model smooth bending. It also allows straight-forward creation of several different levels of detail without having to edit the kinematic joint hierarchy itself.

According to the invention, the default human model combines the kinematic skeleton model with an attached approximation of the spatial extent of the human body model.

More particularly, the approximation may be represented as a sum of un-normalized Gaussians (SoG).

A single Gaussian B has the form:

$\begin{matrix} {{\mathcal{B}(x)} = {\exp\left( {- \frac{{{x - \mu}}^{2}}{2\sigma^{2}}} \right)}} & (2) \end{matrix}$

where σ² is the variance, μ∈

¹ the mean and d=3 for the human model. The Gaussians model only spatial statistics but not any appearance information like color. In order to achieve real-time performance, one need not model the full covariance matrix.

To cover the full 3D model, several spatial Gaussians may be combined into a Sum of Gaussians K:

$\begin{matrix} {{(x)} = {\sum\limits_{i = 1}^{n}\; {\mathcal{B}_{i}(x)}}} & (3) \end{matrix}$

For the human body model, Eq. (3) describes the spatial extent of the human model. In contrast to explicit surfaces, the present model has infinite spatial support, but the influence decreases with the spatial distance to the mean.

An additional color model C=[c_(i)]_(i) may be stored for each sum of Gaussians K, where c_(i) is the color associated with the respective Gaussian B_(i).

The default human body model may comprise 63 Gaussians and 58 joints. Each Gaussian is attached to or associated with a single parent joint of the skeleton, resulting in a default body model K_(m) that is parameterized by the pose parameters Θ of the kinematic skeleton.

This model may be adapted to generate an actor specific body model that roughly represents the shape and color statistics for each person or actor to be tracked.

In order to obtain the pose parameters Θ of the human body model from a given input image, they may be selected such that a two-dimensional projection of the body model becomes essentially equal to a representation of the input image, according to a chosen norm.

The input image may be represented by a sum of Gaussians. Both the image domain Ω∈

² and the 3D human body model may be represented as Sums of un-normalized Gaussians (SoG). In the image domain, Eq. (3) describes the spatial extend of super-pixels that cluster pixels with similar colors as shown in FIG. 3. An additional color model C={c_(i)}_(i) is stored for each SoG K, where c_(i) is the color associated with the respective Gaussian B_(i). 2D-Gaussians.

In order to find an approximation K₁ that represents consistent pixel regions in a given image I, a single Gaussian B_(i) may be created for each image pixel p_(i) and one may assign to each Gaussian the color value c_(i) ∈

³ of the pixel.

In order to avoid an excessive amount of elements and the associated performance penalty, image pixels having similar color may alternatively be clustered into larger regions. The clustering step may use a quad-tree structure for efficiency. Each of the resulting regions may then be approximated using a single Gaussian B_(i).

FIG. 3 shows an input image and a quad-tree structure with average colors used to generate an approximation of the input image by Sums of Gaussians.

For clustering, a threshold ε_(col) may be used to determine which pixels to cluster together (typically set to 0:15). If the standard deviation of colors on a quad-tree node is larger than ε_(col), the node may be subdivided into four sub-nodes, up to a maximum quad-tree depth of typically 8. Each quadratic cluster is then represented by a Gaussian B_(i) where μ is the center of the cluster and σ² is set to be the square of half the side-length of the node. Furthermore, the average color c_(i) of the cluster may be assigned to the Gaussian.

For comparing a 3D body model K_(m) with a 2D image model K_(I), one may define a projection operator Ψ, which projects a given Gaussian {tilde over (B)}_(i):

of the model to a Gaussian B_(i):

²→

. Given a camera C₁ with respective 3×4 camera projection matrix P₁ and focal length f₁, B=ψ({tilde over (B)}) may be defined as the following operations:

$\begin{matrix} {{\mu = \begin{pmatrix} {\left\lbrack {\overset{\sim}{\mu}}^{p} \right\rbrack {x/\left\lbrack {\overset{\sim}{\mu}}^{p} \right\rbrack_{z}}} \\ {\left\lbrack {\overset{\sim}{\mu}}^{p} \right\rbrack_{y}/\left\lbrack {\overset{\sim}{\mu}}^{p} \right\rbrack_{z}} \end{pmatrix}}{s = {\overset{\sim}{s}\; {f_{l}/\left\lbrack {\overset{\sim}{\mu}}^{p} \right\rbrack_{z}}}}} & (4) \end{matrix}$

with {tilde over (μ)}^(p)=P_(l){tilde over (μ)} and [{tilde over (μ)}^(p)]_(x;y;z) being the respective coordinates of the transformed Gaussian mean.

Eq. (4) is only an approximation of the true projection that may be used for computational efficiency. In general, the perspective projection of a sphere is an ellipsoid rather than a circle. However, the error introduced by this approximation proved to be negligible.

According to the invention, a similarity between two given two-dimensional models K_(a):

²→

and K_(b):

²

and the associated color models C_(a) and C_(b), where C={c_(i)}_(i) contains the color values assigned to each Gaussian, may be defined as the integral of the product of K_(a) and K_(b) and a similarity measure between the color models, d(c_(i); c_(j)):

$\begin{matrix} {{{E\left( {_{a},_{b},C_{a},C_{b}} \right)} = {{\int_{\Omega}{\sum\limits_{i \in _{a}}\; {\sum\limits_{j \in _{b}}\; {{d\left( {c_{i},c_{j}} \right)}{\mathcal{B}_{i}(x)}{\mathcal{B}_{j}(x)}\ {x}}}}} = {\sum\limits_{i \in _{a}}\; {\sum\limits_{j \in _{b}}^{\;}\; E_{ij}}}}},\mspace{79mu} {where}} & (5) \\ \begin{matrix} {\mspace{79mu} {E_{ij} = {{d\left( {c_{i},c_{j}} \right)}{\int_{\Omega}{{\mathcal{B}_{i}(x)}{\mathcal{B}_{j}(x)}\ {x}}}}}} \\ {= {{d\left( {c_{i},c_{j}} \right)}2\pi \; \frac{\sigma_{i}^{2}\sigma_{j}^{2}}{\sigma_{i}^{2} + \sigma_{j}^{2}}{\exp\left( {- \; \frac{{{\mu_{i} - \mu_{j}}}^{2}}{\sigma_{i}^{2} + \sigma_{j}^{2}}} \right)}}} \end{matrix} & (6) \end{matrix}$

The similarity d(c_(i); c_(j)) may be modeled by

$\begin{matrix} {{d\left( {c_{i\;,}c_{j}} \right)} = \left\{ \begin{matrix} 0 & {{{{if}\mspace{14mu} {{c_{i} - c_{j}}}} \geq \varepsilon_{sim}},} \\ {{\phi \; 3},{1\left( \frac{{c_{i} - c_{j}}}{\varepsilon_{sim}} \right)}} & {{{{if}\mspace{14mu} {{c_{i} - c_{j}}}} < \varepsilon_{sim}},} \end{matrix} \right.} & (7) \end{matrix}$

with φ_(3;1)(x) being the C2 smooth Wendland radial basis function with φ_(3;1)(0)=1 and φ_(3;1)(1)=0 [H. Wendland. Piecewise polynomial, positive definite and compactly supported radial functions of minimal degree. In: Adv. in Comput. Math., pages 389-396, 1995.]. This functional results in a smooth similarity measure that allows for a certain degree of variation in color matching and is guaranteed to be 0 if the difference between the colors is larger than ε_(sim), which was typically chosen to be 0:15. This particular similarity measure is one-sided, i.e., it only measures the regions that are similar and does not consider or penalize regions that are dissimilar.

When considering the similarity between a projected two-dimensional human body model Ψ(K_(m)) and an image model K_(I), the special properties of the projection operation have to be taken into account. The projection function above ignores possible self-occlusions that may occur when projecting the 3D model onto the 2D image plane. Several Gaussians may be projected onto overlapping 2D positions and thereby contribute several times to the energy function.

FIG. 4 compares a projection with no occlusions (left) with an occlusion situation (right). A top view of the 3D model is shown in the upper row while the lower row shows a view of the 2D model.

Without occlusions, Eq. (5) calculates a correct overlap of a single element. In this example, the color (blue) and the shape are identical, yielding the similarity E_(ii).

If several 3D model Gaussians project to the same screen space coordinate, their contribution is cumulative, yielding a similarity larger than E_(ii), even though two of the model Gaussians should be occluded.

According to the invention, this issue may be resolved implicitly by modifying Eq. (5) to limit the energy that a single image Gaussian can contribute to the total energy of the model:

$\begin{matrix} {{E\left( {_{I},{\Psi \left( _{m} \right)},C_{I},C_{m}} \right)} = {\sum\limits_{i \in _{I}}^{\;}\; {\min\left( {\left( {\sum\limits_{j \in {\Psi {(_{m})}}}^{\;}\; E_{ij}} \right),E_{ii}} \right)}}} & (8) \end{matrix}$

with E_(ii)=πσ_(i) ² being the overlap of an image Gaussian with itself, defining the maximal possible energy it can contribute.

Eq. (8) limits the contribution of a single 2D image Gaussians, yielding the same similarity E_(ii) for both cases.

If an image Gaussian and a projected 3D Gaussian coincide completely the overlap should be maximal. When projecting a second 3D model Gaussian to a nearby location, it cannot contribute more to the already perfect overlap. While this only partially approximates the effects of occlusion, the inventors found that it is sufficient to resolve most ambiguities introduced by occlusions, while analytic derivatives of the similarity function may still be calculated.

In order to estimate the pose-parameters Θ of the kinematic skeleton model from a set of n_(cam) input images I, an energy function E(Θ) may be defined that evaluates how accurately the model described by the parameters Θ represents what may be seen in the images. The function is based on the proposed SoG models and can be used for es-timating the initial actor-specific body model as well as for tracking articulated motion.

The most important part of the energy function according to the invention is measuring the similarity of the body model in the pose defined by Θ with all input images. Given n_(cam) cameras C₁ with respective SoG images (K₁, C₁) and the 3D body (K_(m), C_(m)) parameterized by the pose vector Θ, the similarity function E (Θ) may be defined as:

$\begin{matrix} {{E(\Theta)} = {\frac{1}{n_{cam}}{\sum\limits_{l = 1}^{n_{cam}}\; {E\left( {_{l},{\Psi_{l}\left( {_{m}(\Theta)} \right)},C_{l},C_{m}} \right)}}}} & (9) \end{matrix}$

For the final energy function E(Θ), one may add a skeleton and motion-specific term:

ε(Θ)=E(Θ)−ωw _(l) E _(lim)(MΘ)−w _(a) E _(acc)(Θ)   (10)

E_(lim) (Λ), with Λ=M Θ (1), is a soft constraint on the joint limits and E_(acc) (Θ) is a smoothness term that penalizes high acceleration in parameter space:

${E_{\lim}(\Lambda)} = {\sum\limits_{l \in \Lambda}\; \left\{ {{\begin{matrix} 0 & {{{if}\mspace{14mu} l_{l}^{(l)}} \leq \Lambda^{(l)} \leq l_{h}^{(l)}} \\ {{l_{l}^{(l)} - \Lambda^{(l)}}}^{2} & {{{if}\mspace{14mu} \Lambda^{(l)}} < l_{l}^{(l)}} \\ {{\Lambda^{(l)} - l_{h}^{(l)}}}^{2} & {{{if}\mspace{14mu} \Lambda^{(l)}} > l_{h}^{(l)}} \end{matrix}{E_{acc}\left( \Theta_{t} \right)}} = {\sum\limits_{l \in \Theta_{t}}\; \left( {{\frac{1}{2}\left( {\Theta_{t - 2}^{(l)} + \Theta_{t}^{(l)}} \right)} - \Theta_{t - 1}^{(l)}} \right)^{2}}} \right.}$

where l_(l) and l_(h) are lower and upper joint limits and Θ_(t-1) and Θ_(t-2) the poses of the previous frames.

The weights w_(l) and w_(a) influence the strength of the constraints and were set to w_(l)=1 and w_(a)=0.05 for the majority of the experiments. The impact of the constraints is evaluated in the experimental section.

FIG. 5 shows different actor-specific models (right) derived from example pose images (left).

The default skeleton and body model may be used to estimate an actor specific model in a pre-processing step. A low number of temporally not subsequent, manually segmented multi-view images of example poses may be used as shown in FIG. 5. The four example poses are chosen to articulate a wide range of skeletal joints and therefore allow a relatively accurate estimation of the bone lengths of the skeleton. For each pose represented through a set of multi-view images, the pose parameters Θ are estimated. Additionally, a common set of shape parameters Θ_(shape) may be optimized that define bone lengths as well as the positions and variances of the Gaussian model for a total of 216 degrees of freedom.

The pose parameters Θ may roughly be initialized to correspond to the initial poses manually, and the similarity measure according to Eq. (9) based on the binary color values c_(i) of the silhouette is maximized using a gradient-ascent approach explained below. After the joint optimization, the color images of each pose may be back-projected onto the 3D Gaussian models and the mean color c_(i) for each Gaussian blob may be calculated taking occlusions into account.

Given an image sequence with m frames, one wants to estimate the pose parameters Θ^(t) for each time-step. In each time-step, the parameters Θ^(t) may be initialized by linear extrapolation of the motion in the previous time-steps, i.e., Θ^(t) ₀₌Θ^(t-1)+α(Θ^(t-1)−Θ^(t-2)) with a set to 0.5 for all sequences. The parameters may now be optimized in order to maximize the energy (10). Because of the analytic formulation of similarity measure according to the invention, the analytic gradient ∇ε(Θ) may be calculated efficiently and used it in the optimization procedure.

In order to achieve fast performance, a conditioned gradient ascent may be used for optimizing the energy function. Simple gradient ascent tends to be slow when optimizing energy functions that consist of long narrow valleys in the energy landscape, as it tends to “zig-zag” between opposing walls. To prevent this, it is proposed to introduce a conditioning vector a, into the optimization

Θ_(i+1) ^(t)=Θ_(i) ^(t)+∇ε(Θ_(i) ^(t))∘ σ_(i)   (11)

Here, ∘ is the component-wise Hadamard product of two vectors. The conditioner σ_(i) is updated after every iteration, according to the following rule:

$\begin{matrix} {\sigma_{i + 1}^{(l)} = \left\{ \begin{matrix} {{\sigma_{i}^{(l)}\mu} +} & {{{if}\mspace{14mu} {\nabla{ɛ^{(l)}\left( \Theta_{i}^{t} \right)}}{\nabla{ɛ^{(l)}\left( \Theta_{i - 1}^{t} \right)}}} > 0} \\ {\sigma_{i}^{(l)}\mu^{-}} & {{{if}\mspace{14mu} {\nabla{ɛ^{(l)}\left( \Theta_{i}^{t} \right)}}{\nabla{ɛ^{(l)}\left( \Theta_{i - 1}^{t} \right)}}} \leq 0} \end{matrix} \right.} & (12) \end{matrix}$

Intuitively this conditioning will increase step-size in directions where the gradient sign is constant, and decrease it if the ascent is “zig-zagging”. This proved to reduce the number of iterations necessary to reach a minimum drastically without having to resort to a more complex and more expensive second order optimizer. The inventors chose μ⁺=1.2 and μ⁻=0.5 for all experiments.

At least n_(iter) iterations are performed for each time-step and the iterations are stopped when ∥∇ε(Θ_(i) ^(t))∘ σ_(i)∥<ε, where n_(iter)=10 and ε=0.002 for most experiments.

The inventive method may be implemented using an array of synchronized and calibrated video cameras for obtaining multi-view videos of the actor's movements, and a standard computer. Since the similarity measure is continuous and differentiable, human pose estimation may be solved nearly in real-time even for many camera views. Computation time can further be reduced by parallelizing the evaluation of the double sum in equation 8 using graphics processing units (GPUs) or multi-processor systems.

FIG. 6 shows tracking results of a method according to an embodiment of the invention as skeleton overlay over 12 input images from different camera views.

Experiments of the inventors demonstrate that the inventive method, which does not assume an a-priori background estimate, works even in relatively uncontrolled settings. The motion may be tracked in some challenging situations including characters closely interacting with each other, and occlusion situations, such as a person sitting on a chair at a desk. A quantitative evaluation of the tracking performance demonstrates the reliability of the inventive approach on 17 sequences with over 12000 frames of multiview video. An unoptimized implementation may process a multi-view sequence recorded with 12 cameras at a resolution of 1296×972 pixels at 5-15 frames per second on a standard computer.

More particularly, the sequences were recorded with 12 cameras at a resolution of 1296×972 pixels at 45 frames per second. The quad-tree based image conversion according to the invention was set to a maximal depth of 8 nodes, effectively limiting the used resolution to 162×121 2D Gaussians. This resolution was found to be a suitable compromise between tracking accuracy in all scenes and processing speed. The HSV color space was used for calculating color similarity in all examples.

The sequences were recorded in a room without any special background and cover a wide range of different motions, including simple walking/running, multi-person interactions, fast acrobatic motions, and scenes having strong occlusions by objects in the scene, such as chairs and tables.

The remaining 5 sequences are taken from a different dataset [Y. Liu, C. Stoll, J. Gall, H.-P. Seidel, and C. Theobalt. Markerless motion capture of interacting characters using multiview image segmentation. In CVPR, 2011]. These sequences were recorded with a greenscreen background and show closely interacting characters fighting, dancing, and hugging. Since the inventive approach does not rely on background subtraction, no explicit use is made of the green-screen, but the approach still benefits from the distinct background color. The segmentation is handled implicitly.

To speed up the calculation of the energy function in Eq. (10), all image Gaussians B_(i) may be removed whose color similarity is 0 to all of the body model Gaussians B_(j), as they will not contribute to the energy function at all. Also, all elements that are outside of an enlarged bounding box of the pose in the last estimated frame may be removed.

The inventive method tracked all 17 sequences successfully at 5 to 15 frames per second depending on the scene complexity and the actor's motion speed. To evaluate the range of image resolution one may reliably deal with, and thus how robust the algorithm is to camera resolution, the images of several sequences were downsampled to 81×60 pixels and a single Gaussian was created for each pixel. In this experiment, no significant loss of the tracking quality would be observed.

A standard gradient ascent optimization was compared to the inventive conditioned gradient ascent on the cartwheel sequence shown in FIG. 6. The standard gradient ascent method required on average 65:38 iterations per frame to converge and failed to track the pose of the arms correctly in some frames. The inventive method on the other hand required only 34:56 iterations per frame and tracked the sequence successfully.

Finally, the accuracy of the inventive tracking approach was evaluated and compared to the method of Liu et al. [Y. Liu, C. Stoll, J. Gall, H.-P. Seidel, and C. Theobalt. Markerless motion capture of interacting characters using multiview image segmentation. In CVPR, 2011] on the ground truth multi-person data set kindly provided. The 38 marker positions of the marker system were associated with virtual markers on the inventive body model in the first frame by calculating the closest Gaussian and attaching the marker to its parent joint. Due to the black motion-capture suit and the fast motion, the sequence is difficult to track. The average distance between markers and their corresponding virtual markers on the tracked model of the invention is 44:93 mm with a standard deviation of 27:16 mm. While this is not as accurate as the result reported in [Liu et al.], namely 29:61 mm 25:50 mm, their tracking algorithm employs a laser scanned model and requires several minutes per frame for tracking, while the inventive method tracks the sequence at roughly 6 frames per second.

Without compensating for self-occlusions in the overlap term as in Eq. (8), the sequence fails to track correctly at all, leading to an error of 229:09 mm 252:04 mm. The impact of the weights w_(l) and w_(a) for the joint limits and the smoothness term in Eq. (10) is shown in Tables 1 and 2:

Table 1 shows the effect of the joint limit weight w_(l) on tracking the ground truth dataset. Disregarding this term will lead to physically implausible joint motions, e.g., knees bending backwards.

TABLE 1 Effect of the joint limit weight w₁ on tracking the ground truth dataset. w₁ 0.0 1.0 2.0 error (mm) 48.29 44.93 47.78

Table 2 shows the effect of the smoothness weight w_(a) on tracking the ground truth dataset. Low smoothness leads to jitter in the motion, while high smoothness decreases tracking accuracy.

TABLE 2 Effect of the smoothness weight w_(a) on tracking the ground truth dataset w_(a) 0.0 0.015 0.05 0.1 0.5 error (mm) 46.75 46.33 44.93 46.37 51.74

The present technology has now been described in such full, clear, concise and exact terms as to enable any person skilled in the art to which it pertains, to practice the same. It is to be understood that the foregoing describes preferred embodiments and examples of the present technology and that modifications may be made therein without departing from the spirit or scope of the invention as set forth in the claims. Moreover, it is also understood that the embodiments shown in the drawings, if any, and as described above are merely for illustrative purposes and not intended to limit the scope of the invention. As used in this description, the singular forms “a,” “an,” and “the” include plural reference such as “more than one” or “at least one” unless the context clearly dictates otherwise. Further, all references cited herein are incorporated in their entirety. 

1. A computer-implemented method for tracking an object in a sequence of multi-view input video images, the method comprising the steps of: acquiring a model of the object; and tracking the object in the multi-view input video image sequence, using the model.
 2. The method of claim 1, wherein the model comprises a statistical model representing the shape of the object.
 3. The method of claim 2, wherein the statistical model comprises sums of spatially varying functions chosen such that an overlap or convolution of two or more blobs can be approximated efficiently.
 4. The method of claim 3, wherein the statistical model comprises at least one of radial-symmetric functions, monotone functions, three dimensional Gaussians, hyperspheres, multiquadrics, superquadrics or polyharmonic splines.
 5. The method of claim 2, wherein the model further comprises a kinematic model defining the degrees of freedom of the object.
 6. The method of claim 1, wherein the object is a human body.
 7. The method of claim 1, wherein the model is acquired by scaling, translating or rotating a generic object model.
 8. The method of claim 1, wherein the multi-view input video cameras are synchronized or gen-locked.
 9. The method of claim 1, characterized in that the model comprises attributes characterizing the acquired object.
 10. The method of claim 1, wherein the model is acquired online during tracking.
 11. A system for tracking an object, the system comprising: a range of multi-view video cameras; a memory for storing a model of the object; and a processor for tracking the object based on the model.
 12. A computer-implemented method for acquiring a model of an object, the method comprising the steps of: acquiring a plurality of images of example poses of the object; estimating the model based on the acquired images; and storing the model in a computer-readable medium.
 13. The method of claim 12, wherein the images are multi-view images.
 14. The method of claim 13, wherein the number of multi-view images is low.
 15. The method of claim 13, further comprising the step of manually segmenting the multi-view images.
 16. The method of claim 13, wherein the step of estimating the model further comprises: estimating pose parameters for each example pose of the object; and estimating a set of shape parameters of the object.
 17. The method of claim 16, wherein the pose parameters are estimated using a conditioned gradient ascent.
 18. The method of claim 12, wherein the step of estimating the model further comprises estimating the model based on at least one of marker-based systems, magnetic systems, accelerometer based systems, and motion capture systems.
 19. The method of claim 12, wherein the model is estimated based on an input of at least one of a depth camera, a 2.5D scanner, or a 3D scanner.
 20. The method of claim 12, further comprising the step of manually constraining the model.
 21. A computer-readable medium comprising a model, wherein the model is acquired by performing the steps of: acquiring a plurality of images of example poses of an object; estimating the model based on the acquired images; and storing the model in the computer-readable medium. 